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Diversity-Multiplexing Tradeoff of Double 
Scattering MIMO Channels 

Sheng Yang and Jean-Claude Belfiore 



Abstract 



It is well known that the presence of double scattering degrades the performance of a MIMO channel, 
in terms of both the multiplexing gain and the diversity gain. In this paper, a closed-form expression of 
the diversity-multiplexing tradeoff (DMT) of double scattering MIMO channels is obtained. It is shown 
that, for a channel with rij transmit antennas, tir receive antennas and rig scatterers, the DMT only 
c/2 , depends on the ordered version of the triple (nx, ns, njA, for arbitrary rvj,n% and tir. The condition 

o 

under which the double scattering channel has the same DMT as the single scattering channel is also 

CN| ' established. 
> 



I. Introduction and Problem Description 

Multiple antennas are known as an important means to increase channel capacity and to 
mitigate channel fadings [1], [2]. The tradeoff between the multiplexing gain and the diversity 



CZ3 

■ gain for Rayleigh MIMO channels in the high SNR regime is characterized by the diversity- 



multiplexing tradeoff (DMT) proposed by Zheng and Tse [3]. However, the independent and 



identically distributed (i.i.d.) Gaussian property of the entries of MIMO channels is only estab- 
lished under an idealistic assumption. Recently, a more general model which shows the scattering 
structure in the propagation environment has been proposed [4]. This model considers the rank 
deficiency as well as the fading correlation, by characterizing the channel matrix as a product 
of two statistically independent complex Gaussian matrices. 

The presence of double scattering degrades considerably the performance promised by MIMO 
channels, for both the multiplexing gain and the diversity gain. Intuitively, the performance 
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of double scattering MIMO channels is not better than either the transmitter-scatterers or the 
scatterers-receiver channel. One interesting question is: "what is the impact of double scattering 
on the channel's capability of capturing diversity and providing multiplexing gain in the high 
SNR regime ?". This question is answered in this work, by studying the DMT of double scattering 
MIMO channels. 

More precisely, the main contribution of this work is to provide a closed-form expression 
of the DMT for general double- scattering MIMO channels. It is shown that, for a MEMO 
channel with n T transmit antennas, n$ scatterers and n R receive antennas, the DMT only depends 
on the ordered triple of (n T , ns, n R ). This property can be seen as a generalization of the 
reciprocity of MIMO channels. It is also shown that the upperbound on the channel diversity 
order n T n s n R / max {n T , ns, n R } is usually not achievable, unless for (n T ,n s ,n R ) satisfying 

2 max {n T , n s , n R } + 1 > n T + n s + n R . 

In this paper, we use boldface lower case letters v to denote vectors, boldface capital letters M 
to denote matrices. CM represents the complex Gaussian random variable. [-] T , [-] f respectively 
denote the matrix transposition and conjugated transposition operations. ||-|| is the vector norm. 

means max(0,x). Det(M) is the absolute value of the determinant (let (M). The square 
root P 1 ' 2 of a positive semi-definite matrix P is defined as a positive semi-definite matrix such 
that P = P 1 ^ 2 (P 1 ^ 2 ) t . The dot equal operator = denotes asymptotic equality in the high SNR 
regime, i.e., 

Pi = P2 means hm - — = hm 



snfwoo log SNR snfwoo log SNR 
The rest of the paper is organized as follows. Section II introduces the channel model, 
some preliminaries on complex Wishart matrices and the DMT. Section III studies the DMT 
of Rayleigh product channels, a particular case of the double scattering channel. The DMT of a 
general double scattering channel is provided in Section IV. Section V draws a brief conclusion 
on this work and the Appendix is dedicated to some lemmas and their proofs. 
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II. System Model and Preliminaries 

A. Channel Model 

In this paper, we consider the double scattering MIMO channel with n T transmit antennas, n s 
scatterers and ?t.r receive antennas in the following form 



y = VCSMRHx + z (1) 

with 

h 4 *i /2 fri*J /2 fr 2 *y 2 (2) 

where x E C nT is the transmitted signal with i.i.d. unit variance entries; y E C nR represents the 
received signal; z E C nR is the AWGN with z ~ £A/(0, 1); the constrant C is the normalization 
factor such that SNR is the average Signal to Noise Ratio per receive antenna. Hi E C nRXns 
and H 2 E C nsX " T are statistically independent matrices with i.i.d. unit variance Gaussian entries. 
Correlations at each node are characterized by $ T , <fr s and <& R which are assumed to be positive 
definite matrices 1 with respective dimensions rtj x n T , n s x n s and n R x n R . We denote such a 
channel, a (n T ,ns,n R ) channel hereafter. 

B. Wishart Matrices 

Definition 1 (Wishart Matrix): The m x m random matrix W = HH is a (central) complex 
Wishart matrix with n degrees of freedom and covariance matrix E, (denoted as W ~ W m (n, E)), 
if the columns of the m x n matrix H are zero-mean independent complex Gaussian vectors 
with covariance matrix E. 

Theorem 1 ([5]— [8]): Let W be a central complex Wishart matrix W ~ W m (n, E), where 
the eigenvalues of E are distinct and their ordered values are fi± > . . . > /i m > 0. Let Ai > 
. . . > X g > be the ordered positive eigenvalues of W with q = min{m, n}. The joint p.d.f. of 
A is 

m m , , 

K m , n Det [e~ x ^) J] ^^X^ f[ (3) 
i=i i<j ^ Pi 

'The correlation matrices are positive semi-definite in general. However, it is always possible to have an equivalent channel 
model of positive definite $'s and Gaussian matrices Hi's of reduced dimensions, using the eigenvalue decomposition of the 
correlation matrices and the unitarily invariance property of Gaussian matrices. In this case, the effective numbers of antennas 
and scatterers are rij, n' R and n s , i.e., the respective ranks of $t, 4?r and $s- 
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for n > m, and 



for n < m with 



lib ^ IL 

G m , n Det(E) J] _ HiX, - Xj) 



(4) 



1 ,, ,,m— n— 1 ,,m— n— 1 ~ 777 



1 11 ,,m—n—l ,,m—n—l„~T, — 



^-„-l e - Ml 



, .m—n—1 — - — 



(5) 



i^ m , n and G m ,n are normalization factors. In particular, for E = I, the joint p.d.f. is 



(6) 



i=i 



i<j 



C. Diversity-Multiplexing Tradeoff 

Definition 2 (Multiplexing and diversity gains [3]): A coding scheme {C(SNR)} is said to 

achieve multiplexing gain r and diversity gain d if 

.. P(SNR) logP e (SNR) 
nm ~, — = r anc * hm — - — — — — = —a 
sNR^oologSNR snr^oo logSNR 

where i?(SNR) is the data rate measured by bits per channel use (PCU) and P e (SNR) is the 
average error probability using a maximum likelihood (ML) decoder. 
For any linear fading Gaussian channel 



y = VSNRHx + z 

where z is an AWGN with E{2Z ; } = I and x is subject to the input power constraint Tr {E [xx*]} < 
1, the DMT d(r) can be found as the exponent of the outage probability in the high SNR regime, 
i.e., 

P out (r log SNR) =Prob{logdet (I + SNRHW) <r log SNR} 
= Probjdet (1 + SNRi/JT) < SNR r } 

= SNR- d{r) . (7) 

Lemma 1 (Calculation of diversity-multiplexing tradeoff): Consider a linear fading Gaussian 
channel defined by H for which det (1 + SHRHW) is a function of v, a vector of positive 
random variables. Then, the DMT d(r) of this channel can be calculated as 

d(r) = inf e(a) 

0(a,r) 
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where a { = — log log SNR is the exponent of v„ (D(ot,r) is the outage event set in terms of 
ct and r in the high SNR regime, and is the exponential order of the p.d.f. p a (a) of a, 
i. e. , 

p a (a) = SNR~ £(a) . 

Proof: This lemma is justified by (7) using Laplace's method, as shown in [3]. ■ 
As an example, the DMT of an n R x n T Rayleigh MIMO channel is a piecewise-linear function 
connecting the points (k, d{k)), k — 0,1, ... , min{n R , n T }, where [3] 

d(k) = (n R - k){m - k). (8) 

III. Diversity-Multiplexing Tradeoff of Rayleigh Product Channels 

In this section, we study a special case of the double scattering MIMO channel, where <&t,$s 
and <& R are identity matrices. We call it a Rayleigh product channel. 

Theorem 2: Let H = H 2 H\ with H 2 E C nxl and H 1 E C' xm being independent Gaussian 
matrices with i.i.d. C7V(0, 1) entries. Define (M,N,L) be the ordered version of (m,n,l) with 
M < N < L. Then, the diversity-multiplexing tradeoff of the fading channel 



/SNR 

y=\ -j — Hx + z 

V / m 

is a piecewise-linear function connecting the points (k, d(k)), k = 0, . . . , M, where 



d{k) = (M -k)(N -k) 



[(M-A-k)- 



i2 



(9) 



with A = L - N. 

Before going to the proof, some remarks can be made about the DMT of a Rayleigh product 
channel. 

Remark 1: From (9), we note that 

1) The DMT does not depend on the triple (m, n, I) but only on the ordered triple (M, N, L), 
which can be seen as a generalization of the reciprocity property [2] of MIMO channels; 

2) The DMT of a Rayleigh product channel is always inferior to that of an M x iV Rayleigh 
channel, i.e., d(k) is upperbounded by d(k) = (M — k)(N — k); 

3) The upperbound d(k) is achieved for k > M — A — 1, which means that d(k) coincides 
with d{k) at least for the last section of the curve; 
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4) When L+l > M+N, the Rayleigh product channel has exactly the same DMT performance 
as an M x N Rayleigh channel; 

5) Finally, as a consequence of the previous observation, a Rayleigh product channel is always 
equivalent to an N x 1 Rayleigh channel when M — 1. 

We should point out that the relation between the Gaussian coding bound and the outage bound 
studied in [3] is intimately related to the Rayleigh product channel. In [3], it is shown that the 
Gaussian codeword matrix should be long enough to achieve the DMT of the Rayleigh MEMO 
channel. The code length condition is exactly the same as the condition provided by observation 4 
in the remark above. 

As in [3], the DMT is obtained from the p.d.f. of the eigenvalues ofQ H = HH\ which depends 
on (m, n, I). For now, we know that Q x — H\H\ ~ Wz(m, I). Let us define the eigenvalues of 
Q 1 as iii > . . . > /! m i n {^ m }. Then, Q H = ^H 2 Q 1 / 2 ^j ^H 2 Q 1 i' 2 ^j has the same eigenvalues as 
Q G = GG^ with G = Q\ H\. By definition, conditionned on H\, we have Q G ~ W^r^Qi). 
Therefore, from now on, we can study the eigenvalues A x > . . . > A min {/ imjn }Of Q G , whose joint 
p.d.f. only depends on the eigenvalues of Q x , according to Theorem 1. In the rest of this section, 
we prove Theorem 2 in two cases : min{m, n} > I and min{m, n} < I. 

A. The min{m, n} > I Case 

In this case, we can exchange m and n, by the reciprocity property of MIMO channels. 
Without loss of generality, we assume that m> n. From (3) and (6), we get the joint p.d.f. of 



(A,m) 



i=l i<j 
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where Ci jTri n is the normalization factor. Define a, = — log A*/ log SNR and ft = — log /ij/ log SNR 
for i = 1, . . . , I. Then, we have 

i 

Pa A<*,P) = G, m ,n(logSNR) 2 'nSNR^ (n -' +lK SNR-( m -^ 

i=l 



\ (SNR~ Ql - SNR^)(SNR~ ft - SNR" ft 



i<j 



■ exp (- SNR _/3i j Det exp (-SNR- ( ^- ft) ) 

First, we only consider ft > 0,V2, since otherwise, exp (-EiSNR- ft ) would decay expo- 
nentially with SNR [3]. The high SNR exponent of the quantity Det exp (-SNR~ (Qy " ft) j is 
calculated in Lemma 2. From (20), we only need to consider a.i > ft,Vi, so that ^^(ck, (3) 
does not decay exponentially. Finally, by Lemma 1, the DMT d(r) can be obtained by solving 
the optimization problem 

d(r) = inf e(a,0) (10) 

0(ac,p,r) 



with 



C(a, ft r) 4 | (a, 0) : g(l - + < r, °* ^ " ' ^ ^' , a, > A > 0, V* 
and 

e(a, (3) = - i + + ^(m - n + / - i)ft + ^(a* - A) + - (11) 

i=l i=l i<j 

The optimization problem (10) can be solved in two steps: 1) find optimal (3 by fixing a., and 
then 2) optimize a. Let us start from the feasible region 

< ft = c*i < ft = a 2 < • • • < ft = aj (12) 

in which we have ^^(ai — ft) + = 0. Note that for each j, the feasibility conditions require that 

ft should only move to the left in terms of its positions 2 relative to the ccj's and that A should 
never be on the left of A f° r i > j- Each time A passes an a, from right to left, J2i<j( a i~ 
increases by a, — A, which increases the coefficient of by 1 and decreases the coefficient 
of A by 1. To minimize the value of e(a,/3), A is allowed to pass a { only when the current 

2 The position here refers to the position in the inequality chain of a;'s and /3;'s in increasing order, as the one in (12). 
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coefficient of (3j in (11) is positive 3 . The maximum number of that can be "freed" by j3j is 
j — 1, i.e., ctj-i, Note that the initial coefficient of fij is m — n + I — j and is decreasing 

with j while the number j — 1 is increasing with j. Let j* be the largest number such that 
m — n + l—j > j — 1. Obviously, for j < j*, . . . , «i can be freed and the final coefficients 
of /3j is m — n + l — 2j — 1 (> 0) and /3* = 0. For j > j*, f3j can only free . . . , Oj_( m _ n+ j_j) 
and the final coefficient of /3j is 0. Substituting the optimal solutions /3*'s back into (11), we get 



e a 



1 + Ci)aj 



(13) 



i=i 



where q can be found with the help of Fig. 1(a). Finally, we have 



l—(m—n) 



e la] 



i=l 



( n + 1 - 2i + 



/ + « + (m — n) 



a; 



i 



+ 



E 

i=l— (m— n)+l 



(n + Z + 1 - 2i) cti 



where the coefficient of a» is non-negative and is non-increasing with i. Hence, the optimal 
solution is a* — 1, i = k + 1, . . . , I and a* = 0,i = 1, . . . , k, from which we can verify that 

i2 



d(k) = (l-k)(n-k) 



[(/ - (m - n) - k) + \ 



(14) 



B. The min{m, n} < I Case 

Again, by the reciprocity property, we assume that n < m. However, we should study the 
m > I case and the m < I case separately. We start with the former case. 
1) The n < I < m Case: From (4) and (6), we get the joint p.d.f. of (A, fi) 



P\A\ A*) = B i,m,n Yl^T Yl it* ~ to) 

i=l i<j 



n 



J(A l -A,)Det (S) 



i<j 



3 When the coefficient of /3j in (11) is positive, decreasing f3j decreases e(a,/3). 



(15) 
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g(j) = 2] -Q + m-n) 




2j — (I + m — n) 



(a) The mm{m, n} > / case 



(b) The min{m, n} < I case 



Fig. 1. For each j, the black dots represent the a's that are freed by f3j. For each i, the number of black dots a is the 
coefficient of ou. Thus, for i < g(l), a = {g~ 1 (i)\ - [7 _1 (i)] + 1; and for i > g(l), a = I - |7 -1 (i)] + 1. 



where Bi mn is the normalization factor. Same procedure as the previous case and Lemma 3 
lead to the following asymptotical p.d.f. of (ot, (3) 



Pt 



,/3(a,/3) = (logSNR)' +n ]^[SNR- (ri - l+1)ai 



i=i 



(l+m—n—i)a 



' SNR 



n+l 

JJSNR 

i=l 
n I 

II SNR- (a ^ )+ JJSNR" 

i=l j=n+l «<J 



i=n+2 
n 



(l+m+l-2i)ai 



■ exp ( - SNR~ ft J exp ( - SNR- (a *" ft) 

As before, we only consider $ > 0, Vz, and a» > Pi, for z = 1, . . . , n, in order that p ai p(a., (3) 
does not decay exponentially. Finally, the DMT d(r) can be obtained by solving the optimization 
problem (10) with 

, ai < ■ ■ • < a r 
[L - (my < r, 

i=l 



I 1 ' , > & > 0, for % — 1, . . . , n 
Pi<---< Pi 
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and 



n n+1 

e(a, (3) = y^{ n + 1 ~ i) a i + + m ~ 71 ~ OA 

i=\ i=X 

I n I 

+ £(Z + m + l- 2i)# + ^ ^ («* - 

i=n+2 i=l j'=n+l 



(16) 



The optimization procedure is exactly the same as in the previous case. With the optimal /3/s, 
we have 

n 

e(a) = y^( n — * + 1 + c i) a i 



(17) 



j=i 

where q can be found with the help of Fig. 1(b). Finally, we have 

l—(m—n) 

e(a) = 

i=l 



I + i + (m — n) 



Q; 



+ ^2 {n + l + l-2i)ai 

i=l—(m—n)+l 

where the coefficient of ctj is non-negative and is non-increasing with i. Hence, the optimal 
solution is a* = 1, i = k + 1, . . . , I and ofi = 0, i = 1, . . . , k, from which we have 



d(k) = (l-k)(n-k) 



[(n - (m - Z) - k)+f 



(18) 



2) The n < m < I Case: In this case, pi m +i = • • ■ = Hi = with probability 1. Let 
/*=[//!••■ /i m ] T be the vector of the nonzero eigenvalues of Q P The conditional p.d.f. px\^.(X\fji) 
is given by Lemma 4. The p.d.f. of n being known from (6), we get the joint p.d.f. of (A, fjt) 
in exactly the same form as (15), except that Z and m are interchanged. We have directly 

[(n- (Z -m) - fc)+] 2 



eZ(/c) = (m — k)(n — k) 



(19) 



IV. Diversity-Multiplexing Tradeoff of Double Scattering MIMO Channels 

In this section, we study the DMT of a general double scattering channel, where the antenna 
and scatterer correlations $ T , 4> s and $ R are non-trivial. 

It is intuitive to expect that the DMT is independent of the correlation matrices, as long as 
they are not singular, since the DMT is an asymptotical performance measure. First of all, it is 
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easy to show that the antenna correlations $ T and $ R do not affect the tradeoff. To see this, 
note that 

det(I + SNRHW) = det(I + Sm$ 1 J 2 H 2 $l /2 H l $ J H\$ 1 s /2 H 2 <f>]/ 2 ) 
= det(I + SMR^J 2 H 2 ^l /2 H 1 H\^ 1 s /2 Hl^ 2 ) 
= det(I + SNRHl$l /2 H t 2 H 2 $ 1 s /2 H 1 ) 

where <fr T and <& R disappear in the high SNR analysis. Now, it remains to show that $ s has no 
impact on the high SNR analysis. The following proposition confirms this statement. 

Proposition 1: Let M be any m x n random matrix and T be any m x m non-singular 
matrix whose singular values satisfy cr min (T) = <r max (T) = SNR . Define q = min{m, n} and 
M = TM. Let (Ti(M) > . . . > o q (M) > and a x (M) > ... > cr q (M) > be the distinct 
ordered singular values of M and M, Then, we have 

a,(M) = 0;(M), Vz. 

Proof: For m > n, we consider the left polar decomposition M = UM , where U 
is a m x n matrix with orthonormal columns and M a n x n positive definite matrix with 
<Ti(M) = CTj(M ) for z = 1, ... ,n. Let TU = VT be the left polar decomposition of TU. 
Then, we have o,(M) = ai(T M ) for % — 1, . . . , n. 

For m < n, we make a right polar decomposition M = M W , where U is a n x m matrix 
with orthonormal columns and M a m x m positive definite matrix with cr^M) = o"j(M ) for 
z = 1, . . . , n. Then, we have a^M) = rr ( (T () M () ) for z = 1, . . . , n with T„ = T. 

In both cases, the original problem is equivalent to showing that 

ai(T M ) = ffi(M ), for z = 1, . . . , q, 

with T and M now invertible. Let A and 1? in Lemma 5 be T and M , respectively. By 
applying (28) and (29) of appendix, we have 

a i (M )a m (T ) < <n{M) < a,(M K(T ), 

from which we prove the proposition since a 1 (T ) = a m (T ) = SNR and a^M) = (Tj(M ). ■ 
This proposition says that any invertible transformation with bounded (asymptotically in high 
SNR regime) eigenvalues does not change the asymptotical p.d.f. of the singular values of a 
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random matrix. According to this proposition, we know that the singular values of $ s Hi have 
the same asymptotical p.d.f. as the ones of H±, which leads to the main result of this work. 

Theorem 3: For a (n T , n s , n R ) double scattering MEMO channel (1) with H defined in (2), the 
diversity-multiplexing tradeoff is a piecewise-linear function connecting the points (k, d(k)), k = 
0, . . . , M with d(k) being defined in (9), where (M, N, L) is the ordered version of (n T , n s , n R ) 
with M < N < L. 

Proof: This is a direct consequence of Theorem 2, since the eigenvalues of ($s /2 #i) (Q^hJ 
have the same asymptotical p.d.f. as that of HiH\. ■ 
Note that all observations in Remark 1 apply for the general double scattering MIMO channel. 
In particular, the optimality condition L + 1 > M + N m observation 4 of Remark 1 in terms 

of (nj,n s ,n R ) is 

2 max {n T , n s , n R } + 1 > n T + n s + n R , 

which is also the condition under which the maximum channel diversity order njnsn R / max {n T , n s , n R } 
is achieved. Moreover, this theorem implies that antenna or scatterer correlation does not, indeed, 
have any impact on the DMT of a double scattering channel, as long as the correlation matrices 
are non-singular. Finally, in the singular correlation matrices case, it is straightfoward to show 
that Theorem 3 is still true, but with (riT,ns,n R ) replaced by (n' T ,n' s ,n' R ), the respective ranks 
of the correlation matrices. 

V. Conclusion 

We studied, in this paper, the DMT of a double scattering MIMO channel and showed that, 
as long as the correlation matrices are non singular, it is equal to the DMT of a Rayleigh 
MIMO product channel. This DMT is always lower than the one of a single scattering (n T x n s , 
n s x n R or n T x n R ) MIMO channel and it is equal to that one for certain values of the channel 
parameters. This result is not only interesting for itself, but it also helps to the calculation of 
the DMT of MIMO Amplify-and-Forward [9] cooperative channels as the relayed link can be 
seen as a Rayleigh MIMO product channel. 



February 1, 2008 



DRAFT 



13 



Appendix 



Lemma 2: 



Det [exp f-SNR-to-^l' 

L V / J i,j=i 

= exp ^-^SNR~ (ai " ft) j SNR" 



Proof: 
Let us define Di = Det 



(20) 



D[ = Det 



exp (-SNR- ( ^- ft) 

e -SNFr( Q i-' 3 i>+SNFr (a i-' 3 i> 



and we have 



SNR _(Q i-i"' 3l) +SNR"( Q, !"' 3 i) 



g-SNR-^i-^^+SNR"^!-^' . . . g-SNR^^i-i^^'+SNR-t 01 !-' 3 !' 



e -EiSNR-C a i-^ 



Det 



3 -SNR-( a i-' 3 i ) _ g-SNR-^i"' 9 !) 



g-SNR"^!"' 3 '-!) _ e -SNR-<- a i-P0 



-SNR-^i-i-'' 1 ' _ g-SNR-^i-i"' 3 ') q 



g-SIMR^^-l-ft-i 5 _ g-SNR" < a l-l~^0 q 



= Det 



-SNR-^i-^) 



e -SNR~( a i~' 3 i) / _ g-SNR't"!-' 3 ;) 



-SNR"^-!"^ 5 



g-SNR-^-i-D , J _ e - 



-SNR-("!-' 3 i) 



e -SNR-("i-ft-i) ^ _ e 



-SNR _ ( Q i _ ' 3 i' 



n 



\-(.<*i-Pl) \ ^ 



-SNR-^i-i^i-i' ^ _ e -SNR- (Q i~i" ,3 i ) 



-SNR _ ( <1 i-' 3 i' 



i=l 



where the equations are obtained by iterating the identity SNR a ±SNR 6 = SNR a for a < b. 
Since 1 - e~ x « x for x close to 0+ we have 1 - e -SNR-^-^ = SNR~ K ~ A) if on > fa and 



1-e 



-SNR~( a i -/3 i) 



= SNR otherwise. As shown in the recursive relation above, we must have aij > 



/3i, Vz, in order that A does not decay exponentially. Thus, we have D\ = SNR ^i<^ ai 
and in a recursive manner, we get (20). 
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Lemma 3: 



Det (S) = JJSNR 



n+l I 

(l+m-n-ijoti SNR~^ 



i=l 

n I 



i=n+2 
n 



SNR- (Ql "^ )+ JJSNR' 

i=l j=n+l 



exp -^SNR~ (a! - ft 



Proof: First, we have 



%=i 



Det(S) = JJ^j-^Det 



i=l 



I'l 



(l-n-1) 



■(l-n-1) 



N 

1 e 



-Ai//ii 



-A„/a»i 



(21) 



Then, let us denote the determinant in the right hand side of (21) as .D and we rewrite it as 

■ ■ e~ Ai /^! — e -A iM . . . e --W/ii _ g— a„/mz 



L> = Det 



,(I-n-l) 
tt l,Z 



.(J-n-1) 
a /-lJ 



e" Al/w - 



-Ai/« 



■(I-n-l) 



g— A n /Mi-1 g— An/W 



"An /Mi 



(22) 



Det 



,(I-n-l) 
"l,Z 



(i) 



-Ai//ii 



d 



(Z-n-1) 
i— 1J 



jC 1 ) p -Ai//i(_i 



-An/Ml 



-An/^i 



[(1- e ~ Ai/w ) (23) 



i=i 



where dj^ = ^ B — ^ ft and the product term in (23) is obtained since 1 — e 
1 — e~ Xi ^ 1 for all j < /. Let us denote the determinant in (23) as D[. Then, by multiplying the first 



-(Ai/w-Ai///j) 



column in Di with /xj " 1 and noting that ji\ n i d l -' ; 



l-n-1 j(i-n-l) 



J-n-1 



1, the first column 



of becomes all 1. Now, by eliminating the first / — 2 "l"s of the first column by substracting 
all rows by the last row as in (22) and (23), we have ^j~ n-1 Dj = IliLi ( x _ e~ Al/w ) A-i- By 
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continuing reducing the dimension, we get 

Det(H) = Det[e^/-]^ i n^ n " 1 II £ 



n+l 



n I 



i=l 



i=n+2 



) 



n n (i- 

i=l j=n+l 

from which we prove the lemma, by applying (20). ■ 
Lemma 4: Let W be a central complex Wishart matrix W ~ W m (n,S) with n < m, where 
the ordered eigenvalues of E are jii > . . . > fii > = . . . = /i m = with / > n. The joint 
p.d.f. of the ordered positive eigenvalues Ai > . . . > A n of W equals 



G m , n Det(so n _ - X i) 

i. , !/,; I'll i<j 



(24) 



with 



n A 



1 /ii 



i-n-1 



e m 



I [ii •■■ H t ft e « 



(25) 



Proof: We can prove it by successively applying the l'Hospital rule with /i m , . . . , fj, i+ i — ► 
on the expression in (4). Let us prove by induction that 

Pet (5) DetgQ 

nm n m / r = 3 • (26) 

For I — m — 1, (26) is obviously true. Then, assuming that (26) holds for given I, then, as long 
as I — 1 > n, we have 



lim 



Det(S) 



lim 



DetfS 



Det(5 t _i) 



(27) 



where (27) is deduced from (26). ■ 
Lemma 5: Let A and B be two m x m non-singular matrices. For any n x n matrix M, 
let o"j(M) be the ith largest singular value of M and Tji(M) be the zth smallest one 
(i,(M) = //„ + i_.,(M )). Then, we have 

^■^(AB) < crj(-B) (28) 

rj i+j -i(AB) > Vi(A)vj( B ) (29) 
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for 1 < {i,j} < m and i + j < m + 1. 

Proof: Let AB = UQ be the left polar decomposition of AB with U unitary and Q positive 
definite. Then, we have Q = WAB and Oi(Q) = <7j(AB),Vz. The quadratic form x^Qx can be 
bounded as 

\x f Qx\ 2 = \ (AVx)'iBx)\ 2 

< \\A'Ux\\ 2 \\Bxf (30) 
= (x[Q A xi) (x^Q B x) 

where x\ = Ux, Q A = AA^ and Q B = BB\ The eigenvalue decomposition of Q A and Q B gives 

m 

QA = J2 (J i( A ) ZiZ i and 

i=l 

m 
i=X 

where Zj and are eigenvectors of Q A and Q B , respectively. Now, taking Xk = U^Zk for 
k — 1, . . . , i— 1 andxfc = for A; = z, . . . , i+j—2, we have, Vx J_ Xk for = 1, . . . , i+j—2, 

iUx) Q A (Ux) < aj(A) \\x\\ 2 (31) 
x*Q B x<a?(B)\\x\\ 2 . (32) 
From (30), (31) and (32) and the Courant-Fischer theorem [10], we get 

o*(AB)< max < a 2 (A)a 2 (B), 

XXX-L,...,X i+ j-2 \\x\\ 

from which we have (28). 

Note that for any invertible matrix M, we have rji(M) = a~ 1 (M~ 1 ). By applying this equality 
and using the inequality (28), it is straightfoward to get (29) after some simple manipulations. 
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